Empirical traffic data and their implications for traffic modeling 



Dirk Helbing 

//. Institute of Theoretical Physics, University of Stuttgart, Pfajjenwaldring 57/III, 70550 

Stuttgart, Germany 

Abstract 

From single vehicle data a number of new empirical results about the tem- 
poral evolution, correlation, and density-dependence of macroscopic traffic 
quantities have been determined. These have relevant implications for traffic 
modeling and allow to test existing traffic models. 
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With the aim of optimizing traffic flow and improving today's traffic situation, several 
models for freeway traffic have been proposed, microscopic and macroscopic ones @-|14 



Only some of them have been systematically derived from the underlying laws of individual 
vehicle dynamics |Tl|-p!^]. Most models are phenomenological in nature 0-0. These base on 



various assumptions, the correctness of which has not been carefully discussed up to now, 
mainly due to a lack of empirical data or difficulties to obtain them. Therefore, this paper 
presents some fundamental empirical observations which allow to test some of the models. 

The empirical relations have been evaluated from single vehicle data of both lanes of 
the Dutch freeway A9 between Haarlem and Amsterdam (cf. Fig. [l|) These data 

were detected by induction loops at discrete places x below the lanes i of the roadway and 
include the passage times ta{x, i), velocities Va{x, i) and lengths i) of the single vehicles 
a. Consequently, it was possible to calculate the number Ni{x, t) of vehicles on lane i which 
passed the cross section at place x during a time interval [t, t + T], the traffic flow 

Q,{x,t) ■.= N,{x,t)/T, (1) 

and the macroscopic velocity moments 

(-')^--=WI7T^ E K{x,^)r. (2) 

If nothing else is mentioned, the interval length was chosen T = 5min, since this allowed 
to separate the systematic temporal evolution of the macroscopic traffic quantities from 
their statistical fluctuations |jT6|. The vehicle densities pi{x,t) were calculated via the flow 
formula 

Qi{x,t) = pi{x,t)Vi{x,t) . (3) 
A detailled comparison with other available methods for the determination of the average 



velocities Vi := {v)i and the vehicle densities pi from single vehicle data will be given in |]T6 
Finally, the lane averages of the above quantities were deflned according to 
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J2 E Ki^,^)]\ 



(5) 



i t<ta<t+T 



and 



p{x,t) := Q{x,t)/V{x,t), 



(6) 



where / denotes the number of lanes and V{x,t) := (v). Therefore, we have the following 
relation: 



For reasons of simplicity, most macroscopic traffic models describe the dynamics of the 
total cross section of the road in an overall manner by equations for the density p and the 
average velocity V. However, one would expect that a realistic description requires a model 
of the traffic dynamics on the single lanes and their mutual coupling due to overtaking and 



oscillations among the lanes p8[| . 

In order to check this, we will investigate the correlation between neighboring lanes. 
Figure ^ shows that the temporal course of the densities pi{x,t) and p2{x,t) is almost 
parallel. A similar thing holds for the average velocities Vi{x,t) and V2{x,t) [|T^. The 
difference between the curves is mainly a function of density (cf. Figure |^): At small 
densities, the vehicles can move faster on the left lane than on the right one, whereas at 
high densities the left lane is more crowded than the right one. In addition. Figure § shows 
that the variances 



behave almost identically on the neighboring lanes (although the order of magnitude of the 
average velocities Vi{x,t) changes considerably). This strong correlation between neigh- 
boring lanes probably arises from overtaking and lane-changing maneuvers. It justifies the 
common practice to describe the dynamics of the total cross section of the road in an overall 
manner. 




(7) 



lane-chaning maneuvers p!3|JT6|] . This could cause a more complex dynamics like density 




(8) 
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Now we face the question, how a reahstic traffic model must look like. Due to the 
conservation of the number of vehicles, the dynamics of the vehicle density is given by the 
continuity equation [P-|6| jr3|JT6[ 



^ + ^ = .+ (..t)-.-(..(), (9) 

where z/+(x, t) and z/^(x, t) are the rates of vehicles which enter or leave the freeway at on- 
and off-ramps, respectively. Lighthill, Whitham and Richards have suggested to specify the 
flow Q{x,t) in accordance with an empirical flow-density relation Qe{p) 

Q{x,t) = QeipM) . (10) 

This relation has been called into question, since the resulting model cannot describe the 



emergence of phantom traffic jams or stop-and-go traffic [Ty|T^. Therefore, some researchers 
have introduced an additional dynamical equation for the average velocity V{x,t) which al- 
lows to describe instabilities of traffic flow 0-|IO|,|14[. However, others have interpreted these 
phenomena as effects of fluctuations or of phantom bottlenecks caused by slow, overtaking 
vehicles like trucks Hence we check relation (|r^) in Figure |^. It is found that (|TUp be- 



comes invalid above a density of about 12 vehicles per kilometer and lane, where a hysteresis 
effect occurs I^TJ. This indicates a transition from stable to unstable traffic flow. 



An empirical proof of emerging stop-and-go traffic is presented in Figure ^. During the 
rush hours between 7:30 am and 9:30 am, average velocity breaks down at place x = 41.75 km 
because of the on-ramp at x = 41.3 km. Nevertheless, the traffic situation recovers at the 
successive cross sections, i.e. average velocity increases again. In spite of this, the initially 
small velocity oscillations at x = 41.75 km grow considerably in the course of the road. This 
corresponds to emerging stop-and-go traffic (i.e. alternating periods of acceleration and 
deceleration). At the same time, the wavelength of the oscillation increases. This is in good 
agreement with computer simulations which show a merging of density clusters leading to 
larger wave lengths |^. 

After we have found that we need a dynamic velocity equation for an adequate description 
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of the spatio-temporal evolution of traffic flow, we have to clear up the question, whether 
we also need a dynamic equation for the variance 

e:={{v-vr) = e + {{v,-vr) (11) 

or not. Theoretical considerations on the basis of gas-kinetic approaches have shown that 
the velocity equation depends on the variance, for which a separate equation can be derived 
||TTHT6|. Nevertheless we will try out the equilibrium approximation 

e(x,t) = ee(p(x,t)), (12) 

where Oe(p) is the empirical variance-density relation (cf. Figure 0). Figure ^ shows that 
this approximation fits the temporal evolution of the variance in a satisfactory way as long 
as the average velocity V does not rapidly change. However, when the velocity breaks down 
or increases, the variance shows mysterious peaks. These are a consequence of having built 
the temporal averages 



{V''){x,t) ^ {V'^)a{x,t) := -J^ dt' {V')a{x,t') (13) 

over finite time intervals T, where {v^)a{x,t) := / dvv'^ Pa{v; x , t) with the actual velocity 
distribution Pa{v; x,t). In linear Taylor approximation we find 

v(x, t) « [v;(i, i) + ^Yi^it' - 1) 

= V.M) + |^. (14) 

Since dVa/dt is varying around zero, the measured value V fiuctuates around the actual 
value Va{x,t) := {v)a{x,t). For the variance we find 



Q={{V- V)^)a = Qa + [Va - V]^ 

Therefore, time averaging leads to a positive correction term which becomes particulary 
large, where the average velocity changes rapidly, but vanishes in the limit T — > 0. This 
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correction term describes the variance peaks in Figure ^ quite well. Consequently, the 
dynamics of the variance can be reconstructed from the dynamics of the vehicle density 
p{x,t) and the average velocity V{x,t). 

Summarizing our results, we were able to demonstrate the following by empirical data: 
1. The dynamics of neighboring lanes is strongly correlated so that the total freeway cross 
section can be described in an overall way. 2. There is a transition from stable to unstable 
traffic flow at a critical density pcr of about 12 vehicles per kilometer and lane. 3. Emergent 
stop-and-go traffic exists, so that a realistic traffic model must contain a dynamic velocity 
equation. 4. The variance can be well approximated by an equilibrium relation, if corrections 
due to time averaging are taken into account. These conclusions seem to be also valid for 
other stretches of freeway systems, at least European ones. 

The empirical findings question the fluid-dynamic model by Lighthill, Whitham and 
Richards 0,^. They are in favour of the phenomenological models by Payne 0, Phillips 



11|, Kiihne m, Kerner and Konhauser M, Hilliges M as well as a recent model by Helbing 



l^ , p!6[] which has been systematically derived from the microscopic vehicle dynamics via a 
gas-kinetic level of description. The last of these models flts the instability region best, in 
particular the surprisingly low critical density Pcr [RE • 
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FIG. 1. The investigated stretch of the Dutch two-lane freeway A9 from Haarlem to Amsterdam 
including on- and off-ramps. Detectors are indicated by vertical lines. The detector at x = 41.3 km 

( ) only evaluates on-ramp traffic and a bus lane. Between x = 40.8 km and x = 37.6 km traffic 

flow is not disturbed over more than three kilometers. The speed limit is 120km/h. 
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FIG. 2. The temporal course of the average velocities Vi ( — : right lane; : left lane) and 

the vehicle densities pi ( — : right lane; • • •: left lane) on October 14, 1994 at x = 41.75km. 
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FIG. 3. Average differences between the vehicle densities (O) and the average velocities (+) on 
both lanes of the Dutch freeway A9 on November 2, 1994 (T = 1 min). 
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FIG. 4. The temporal course of the average velocities Vi ( — : right lane; : left lane) and the 

standard deviations -y/Ol of vehicle velocities ( : right lane; • • •: left lane) on November 2, 1994 

at X = 41.75 km. 
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FIG. 5. Comparison of the fundamental diagram Qe{p) (i-e. the average flow-density relation, 
-) with the temporal evolution of traffic flow Q{x,t) = p{x,t)V{x,t) ( — ) at x = 41.75 km. 
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FIG. 6. Temporal evolution of the average velocity V(x, t) at successive cross sections of the 

freeway on October 14, 1994 (•••: x = 41.75km; : x = 40.8km; : x = 39.6km; — : 

X = 37.6 km). 
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FIG. 7. Average density-dependence of the standard deviation of vehicle velocities (for 
T = Imin, O) and corresponding fit function ( — ). 
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FIG. 8. Temporal evolution of the standard deviation y/Q{x, t) of vehicle velocities ( — ) in 
comparison with the equilibrium approximation ^/Qe{p{x,t)) (•••)• Obviously, the empirical vari- 
ance Q{x,t) has peaks where the average velocity V{x,t) ( ) changes considerably. In order to 

describe these, a correction term due to time averaging must be added ( — ). (Data from November 
2, 1994 at X = 41.75 km) 



